Stable ferromagnetism and doping induced half-metallicity in asymmetric graphene nanoribbons 
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We propose a class of graphene nanoribbons showing strong intrinsic ferromagnetic behavior due to their 
asymmetry. Such ribbons are based on a zig-zag edged backbone surmounted by a periodic, triangular notched 
region of variable size. The electronic properties as a function of the topology are investigated. Interestingly, 
substitutional doping by boron or nitrogen induces half-metallicity. The most effective doping sites can be 
inferred from the band structure. Given the present rapid development of bottom-up strategies for the synthesis 
of atomically precise carbon nanostructures the proposed class of nanoribbons emerges as a real candidate for 
spintronic applications at ambient temperature. 

PACS numbers: 61.48.Gh, 73.22.Pr , 81.05.ue 



Graphene, with its outstanding properties d], could play an 
important role for spin- sensitive devices given the long life- 
time |2] and coherence length ^ of spin-injected electrons. 
Nevertheless, graphene is a semimetal, i.e. a zero bandgap 
material with negligible density of states at the Fermi energy, 
which limits its suitability for electronic or optoelectronic de- 
vices. This can be overcome by structuring graphene at a very 
small size where quantum confinement and edge effects can 
be used to design appropriate electronic band gaps. In partic- 
ular graphene nanoribbons (GNR) are being investigated and 
considered as promising candidates . The recently devel- 
oped bottom-up synthesis methods [5] based on appropriately 
functionalized molecular precursors meet the requirement of 
precision at the atomistic scale and motivate the theoretical 
research on design of GNRs with desired properties. 
The properties of GNRs are mainly determined by the edge 
shape and GNR- width [4]. Unpaired electron states would 
occur at the edges as revealed by Klein et al. |0] based on 
considerations from resonance theory. There are two types of 
edges, armchair and zigzag. Ribbons with zigzag edges show 
localized edge states with energies close to the Fermi level 
10]. These are responsible for strong magnetic response. For 
smooth zigzag GNRs the net spin density on one edge is usu- 
ally balanced by a spin density of the same amplitude on the 
other, and no effective magnetic moment results. In ribbons 
with armchair edges, on the other hand, edge states are absent 
and no magnetism is found 180 . 

Nevertheless, zigzag GNRs could be useful as spin filters in 
spintronic devices due to higher magneto-resistance and 
to a spin-correlation length of at least 1 nm. Realistic den- 
sity functional theory calculations [10] classify zigzag GNR 
as "narrow gap semiconductors" with similar gaps for up 
and down spins. The application of a large external electric 
field makes this system half-metallic, but generation of a field 
larger than 0.1 eV/A could raise stability issues in the system. 
Many efforts have been made to design anti-ferromagnetic 
GNRs. Unconventional methods like rolling [11] and inser- 
tion of line defects 11211 were examined. However, in these 
cases magnetism is sensitive to the the ribbon width and ge- 
ometry. Theoretical studies were done for a number of edge- 




FIG. 1. The (n z = 2; p T = 3) TA-GNR with the A and B sublattices 
colored in green and purple, respectively. Also indicated by A, A , 
A" and B are the doping sites discussed in the text. The colored 
hexagons define the size of the zigzag backbone and of the triangular 
section. 



functionalizations (N,NH 2 , N0 2 , CH 3 ) of zigzag GNR (see 
e.g. Ref. [13]). Other efforts were invested in doping e.g. 
by nitrogen and boron HO] which can induce half-metallicity 
or by Al turning the material magnetic fl 1 411 . Moreover, Yu et 
al. H 1 311 underline the crucial role of the doping position. Anti- 
ferromagnetism is found for the ground state for quantum dots 
with zigzag edg es 1 161 and for hydrogenated derivatives of the 
same geometry II 1 Vfl . 

By means of combined tight binding (TB) and density func- 
tional theory (DFT) calculations, we propose in this paper 
a class of strongly asymmetric GNRs showing stable ferro- 
magnetic behavior as a consequence of their topology. It is 
shown that substitutional doping of the GNRs at particular 
sites induces half-metallicity in the structure. Based on the 
electronic band structure of the undoped GNRs, we finally 
explain why the doping at such sites is efficient with respect 
to half-metallicity. The stability of the ferromagnetic states 
is demonstrated by computing the cost of creating an anti- 
ferromagnetic boundary. 

The main concept is the addition of a triangular- shaped 
structure on one side of a zigzag GNR. The result is a 
triangular asymmetric GNR, dubbed TA-GNR made up by 



2 



a "zigzag backbone" whose width is defined by the number 
nz of benzene rings in the vertical direction, surmounted by 
a triangular shaped section. Both edges are mainly zigzag. 
We include no space between the triangular-regions. The 
size of the triangular section is defined by the number pr 
of phenyl rings on its side and we characterize our GNRs 
through the pairs (nz, Pt)- For example, the ribbon in Fig. [T] 
has a nz - 2 backbone and p? = 3 phenyl rings on the side of 
the triangular section. The proposed design allows to enhance 
magnetism in the ribbon: on the one hand the zig-zag edge 
aims for localized edge states, whereas the triangular section 
introduces an asymmetry between spin-up and spin-down 
density leading to a ferromagnetic ground state. 
The quantum-Espresso code is used to perfor m sp in-polarized 
density functional theory (DFT) calculations 111 811 W19I1 includ- 
ing full optimization of atomic positions and unit cell. For all 
topologies, we observe considerable corrugation mainly due 
to the steric repulsion between hydrogens at the boundary 
point of two triangular sections. For each GNR a series of 
calculations is performed with a constrained magnetization 
between and 10 pb- The resulting wave functions are then 
relaxed completely and the ground state is identified by a 
comparison of the energies. The ferromagnetic ground state 
is reliably identified by this procedure as shown by the large 
stabilization energy of that state (see Table [J). 
Single band, nearest neighbor tight binding (TB) calculations 
with a mean field approximation to the Hubbard Hamiltonian 
lEoll are used to screen a wide range of geometries and 
magnetizations and to better investigate the character of the 
bands for different geometries. 

At the TB level, the band structure of the TA-GNR is charac- 
terized by a large gap (E gap = t for nz = 1 , being t the nearest 
neighbor hopping parameter, t & 3 eV for GNRs) formed by 
a series of weakly spin- split bands. Within this gap we find 
a number of fully spin polarized bands near the Fermi level 
(E F = 0) (Fig. [2). Only the bands of one spin orientation are 
occupied turning the TA-GNRs magnetic. Upon increase of 
the size of the backbone, the gap-width decreases (the zigzag 
metallic nature starts to prevail) and the fully polarized bands 
are modified. Indeed, some of them of them remain flat (three 
for p T = 5 as in Fig. [2]), whereas the remaining bands (two 
for p T = 5) start to show considerable dispersion toward the 
limit of infinite width where valence and conduction bands 
touch at E F 

Fig. [3] shows the almost flat bands in detail for the case 
(nz = 6,Pt = 5). The lowest three bands correspond to 
localized edge states rapidly decaying toward the TA-GNR 
bulk and are scarcely influenced by the ribbon width. The 
upper two bands in turn (states 1, 2, 6 and 7) are delocalized 
and can be considered as TA-GNR bulk bands. A backbone 
makes them disperse (as e.g for nz=34 in Fig. [2]). The plots 
of the wave function indicate that the shifted bands are of 
delocalized character, with considerable amplitude in the 
backbone region (see states 1 and 2). The bands hosting states 
3, 4 and 5 on the other hand have states at the T-point which 
are localized at the zigzag-edges. Departing from the T-point 



the localized character of the same band persists (see states 8, 
9 and 10). At the zone-boundary even the bands hosting states 
6 and 7 become more localized. We have also investigated the 
effects of an increase of the size of the triangular section. We 
compute the net magnetization for the backbone nz = 6 and 
shapes ranging from p T = 1 to pt = 9 at the TB level. We 
find p = IpB to p = 9 pb per unit cell, i.e., a direct correlation 
between p T and the magnetization. 

Although TB calculations can give reliable information 
on the spin-resolved band structure they can not answer if 
the ferromagnetic state is the ground state. To identify the 
energetically lowest lying state we perform DFT calculations 
for selected combination of nz and pt- We compute the 
total magnetization of the most stable structure, together 
with the energy difference between the magnetized and 
non-magnetized state, interpreted as magnetic stabilization 
energy (Table U). In general the stabilization is stronger 
for nz = 1 than for nz = 2. The stabilization energy per 
Pt (note: pt is proportional to the length of the unit cell) 
is approximately constant for nz = 1 (~70 meV) and for 
nz - 2 (^45 meV), except for the narrowest ribbons. The 
magnetization increases with p T from p = 0.97 pb for p T = 1 
to p = 4.86 ps for p T = 5 and nz = 1 (TB value: ptb = 5 Pb)- 

The band structure for (nz = 2;p z = 3) from DFT is given 
in Fig. |4(a)[ It confirms the features found by TB in undoped 
TA-GNRs with nz = 1, 2: The ribbons appear to be insulators 
with a gap of around 1 eV. The flat bands are approximately 
symmetric with respect to the Fermi-level. The number of 
fully spin polarized and at the same time occupied bands 
corresponds to the "amplitude" pr of the asymmetry. The 
further bands split by a gap of around 3 eV. 

TB and DFT agree and highlight an anti-ferromagnetic 
nearest neighbor coupling, turning the lattice bipartite. In the 
following we will refer to the network of carbon atoms starting 
from hydrogen saturated carbons as sublattice A, the remain- 
ders as sublattice B (Fig. Q]) Each sublattice has a preferred 
spin orientation. In all cases the flat bands at the Fermi level 
(Fig. |4(a)| ) consist largely of a superposition of 7r-electrons of 
sublattice A. The confirmation of the magnetic character com- 
ing from DFT simulation leads us to another kind of symmetry 
breaking, namely asymmetry between spin-up and spin-down 
bands induced by doping. 

We substitutionally dope the GNRs with nitrogen and 
boron. The doping removes one electron from the valence 
band for boron, and occupies an additional band for nitrogen. 
We study the effect of substitution at sites A, B, A' and A" 
(see Fig. [T]): the substitution at site A' shifts one flat band. 
For nitrogen this is from the unoccupied upper to the occu- 
pied lower set of three bands. The gap and the flatness of the 
bands remain unchanged. The substitution at site B does not 
affect electronic orbitals contributing to the bands near E F . It 
only shifts the Fermi level into the localized bands leaving the 
bands near E F flat. Although half-metallic, the material is ex- 
pected to be a bad conductor. The substitution at site A shifts 
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FIG. 2. Series of TB band structures in the first Brillouin zone for 
Pt = 5 with n z = 2, 6 and 34. Blue marks spin-up and red spin down 
bands. The magnetization is 5 p B per unit cell in all cases. The on- 
site Coulomb repulsion parameter U equal to the nearest-neighbor 
hopping t has been used. 
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TABLE I. The lattice vector ao, the total magnetic moment, the sum 
of the absolute values of the magnetic moments and the energy differ- 
ence AE = (E - E p )lp T between the unpolarized and the polarized 
state from DFT calculations. 



the Fermi-level into the set of flat bands and at the same time 
induces a considerable dispersion in two bands. The effect on 
the lower band is direct, through re-hybridization of the donor 
(N, Fig. |4(b)| ) or acceptor (B) state. This re-hybridization was 
quantified in Ref. lElll for impurities in symmetric GNRs. The 
effect on the corresponding band in the other spin channel is 
indirect and communicated through the Hubbard repulsion of 
the localized orbitals these bands correspond to. The material 
can turn into a half-metal if one part of the dispersing band 
remains close to the bands of the same spin as it is the case for 
the nitrogen-doping. The effects on the band structure through 
doping at site A" are similar with those of site A, the latter case 
is illustrated in Fig. |4(b)| 

We conclude that the most efficient substitutions are at sites 
A and A" close to either edge, where the density of the edge 
state is not attenuated too strongly. The doping introduces a 
large dispersion while leaving one part of the dispersing band 
unchanged. The resulting material is a half-metal. The reason 
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FIG. 3. Detailed TB insight into the five frontier states of the TA- 
GNR with (n z = 6,p T = 5). The red lobes depict the positive phase 
whereas the blue lobes depict the negative phase. The size of the 
cirles is proportional to the occupancy of that site. The unit cell is 
highlited by a blue rectangle. In the band structure the unoccupied 
counterpart - which would result by mirroring the occupied part at the 
Fermi-level, a result of the particle-hole symmetry — is not shown. 



for the dispersion is following: In general the bands around 
E F localized mostly on the smooth ribbon edge at the T-point 
extend into the central part of the ribbon at the edge of the 
Brillouin-zone boundary (see also Fig.0 states 5 and 10 in TB 
and Fig. |4(a)| ). For undoped ribbons this introduces no disper- 
sion because the localized state feels the same potential in the 
whole Brillouin-zone. In doped ribbons, the states at the zone 
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boundary overlap more strongly with the dopand than those 
in the center. In the case of nitrogen-doping the potential is 
attractive (N +6 ) and the bands are bent downwards at the zone 
boundary (Fig. |4(b)| ). For B doping (Fig |4(c)| ) the bending is 
opposite and bands are also affected by a rehybridization of 
the edge states. This causes an overlap with the dopand for 
the whole Brillouin-zone. This results in a additional rigid 
shift of the band. The bands near E F are only affected by 
substitutions in the A-sublattice. The states around E F have 
no density on the B-sublattice and therefore their interaction 
vanishes with impurities on the B-sublattice. For impurities 
placed in the edge no re-hybridization takes place as the edge 
states are typically located at the impurity site. 

Coming back to the stability issue, we predict the nanorib- 
bons with low values of nz to be ferromagnetic at room tem- 
perature. The difference of the ground state energy Eg of a 
supercell and the energy En of the same cell with half of the 



spins on the A sub-lattice flipped [22] leads to an estimate of 
the nearest-neighbor coupling term / in a lD-Ising model. We 
find J = 0.127 eV (for n z = 1) and / = 0.99 eV (n z = 2), lead- 
ing, in the mean-field approximation, to a Curie temperature 
of T c - 1000 K (n z = 1) and 700 K (n z = 2), well above room 
temperature ll23ll . 

In conclusion, we have introduced a novel class of asymmet- 
ric nanoribbons that show a stable magnetization, and, upon 
appropriate doping, exhibit half-metallic behavior. The lat- 
ter is induced whenever the doping is close to the edge of 
the ribbon on the same sublattice as the hydrogen- satured 
carbon atoms. The band structure is related to the geome- 
try, and the half-metallic properties are potentially technologi- 
cally attractive particularly for spintronics applications: ferro- 
magnetic atomic coupling lEill . magnetic domains and Curie 
temperature above room temperature lEill . As previous re- 
search by our group has shown, novel bottom-up strategies 
are rapidly emerging. They allow the synthesis of graphene- 
derived structures with atomic precision ^ and announce the 
possibility of doping the GNR at a well-defined position. The 
present work thus acquires technological relevance within a 
strategy for the realization of nanomaterials with designed 
composition, size, shape and properties which in turn could 
lead to applications in quantum computation and spin manip- 
ulation. 
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FIG. 4. The band structure and the density of state of the (n z = 
\;p T = 3) ribbon without doping |(a)[ compared to the one with ni- 
trogen [(b)] and boron doping [(c)] The \if/ 2 \ of the state leading to 
half-metallicity is given at different points in reciprocal (isosurface 
at 0.0025 e/cfy. The band leading to half-metallicity changes its 
character across the Brillouin zone. The effect of the dopand on this 
band is more pronounced in the region of the Brillouin zone where 
the overlap of the corresponding state and the dopand is large. |(a)|in- 
set: the magnetic ordering in TA-GNRs (blue: positive polarization, 
red: negative polarization, white: neutral). 



